
R : Copyright 2003, The R Foundation for Statistical Computing
Version 1.8.1  (2003-11-21), ISBN 3-900051-00-3

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for a HTML browser interface to help.
Type 'q()' to quit R.

[Previously saved workspace restored]

> invisible(options(echo = TRUE))
> # estimate John Jackson's model using multinomRob
> #
> 
> options(width=130);
> 
> library(multinomRob)
Loading required package: rgenoud 
Loading required package: MASS 
Loading required package: mvtnorm 
> 
> # read John Jackson's replication dataset (PA 2002, Winter)
> #
> jjdat <- read.table(file="pareplic.asc", header=T);
> 
> #> names(jjdat)
> # [1] "Voivodship"    "ID"            "UD.KLD"        "SLD"
> # [5] "PSL"           "UP"            "Catholic"      "Other"
> # [9] "Constant"      "New.Jobs"      "Unemp"         "SOE"
> #[13] "Church.Att"    "Schooling"     "Age"           "Farmer"
> #[17] "Village"       "Suchoka.UD"    "Pawlak.PSL"    "Kwasn.SLD"
> #[21] "Zioqkowska.UP" "Voters"
> 
> jjuse <- jjdat;
> #total number of voters
> jjuse$TotalVotes   <- jjdat$Voters * 1000;
> jjuse$SLD      <- jjdat$SLD/100      *jjuse$TotalVotes;
> jjuse$PSL      <- jjdat$PSL/100      *jjuse$TotalVotes;
> jjuse$UP       <- jjdat$UP/100       *jjuse$TotalVotes;
> jjuse$Catholic <- jjdat$Catholic/100 *jjuse$TotalVotes;
> jjuse$Other    <- jjdat$Other/100    *jjuse$TotalVotes;
> jjuse$UD.KLD   <- jjdat$UD.KLD/100   *jjuse$TotalVotes;
> 
> jjuse$Age  <-  jjdat$Age /10;
> 
> ## Set parameters for Genoud
> zz.genoud.parms <- list( pop.size             = 10000,
+                         wait.generations      = 30,
+                         max.generations       = 100
+                         )
> 
> # estimate the model as done in the Mebane-Sekhon AJPS 2004 paper
> 
> muljj <-
+   multinomRob(
+     list(SLD ~      New.Jobs + Unemp + SOE + Church.Att + Schooling + Age + Kwasn.SLD ,
+          PSL ~      New.Jobs + Unemp + SOE + Church.Att + Schooling + Age + Farmer + Village + Pawlak.PSL ,
+          UP ~       New.Jobs + Unemp + SOE + Church.Att + Schooling + Age + Zioqkowska.UP ,
+          Catholic ~ New.Jobs + Unemp + SOE + Church.Att + Schooling + Age ,
+          Other ~    New.Jobs + Unemp + SOE + Church.Att + Schooling + Age ,
+          UD.KLD ~   0 + Suchoka.UD ),
+     jjuse,
+     multinom.t = 0,
+ #    genoud.parms = zz.genoud.parms,  # use the defaults
+     print.level = 3);

Your Model (xvec):
                                                                       SLD PSL UP Catholic Other UD.KLD
(Intercept)/(Intercept)/(Intercept)/(Intercept)/(Intercept)/Suchoka.UD   1   1  1        1     1      1
New.Jobs/New.Jobs/New.Jobs/New.Jobs/New.Jobs/NA                          1   1  1        1     1      0
Unemp/Unemp/Unemp/Unemp/Unemp/NA                                         1   1  1        1     1      0
SOE/SOE/SOE/SOE/SOE/NA                                                   1   1  1        1     1      0
Church.Att/Church.Att/Church.Att/Church.Att/Church.Att/NA                1   1  1        1     1      0
Schooling/Schooling/Schooling/Schooling/Schooling/NA                     1   1  1        1     1      0
Age/Age/Age/Age/Age/NA                                                   1   1  1        1     1      0
Kwasn.SLD/Farmer/Zioqkowska.UP/NA/NA/NA                                  1   1  1        0     0      0
NA/Village/NA/NA/NA/NA                                                   0   1  0        0     0      0
NA/Pawlak.PSL/NA/NA/NA/NA                                                0   1  0        0     0      0

[1] "multinomRob:  WARNING.  Limited regressor variation..."
[1] "WARNING.  ... A regressor has a distinct value for only one observation."
[1] "WARNING.  ... I'm using a modified estimation algorithm (i.e., preventing LQD"
[1] "WARNING.  ... from modifying starting values for the affected parameters)."
[1] "WARNING.  ... Affected parameters are TRUE in the following table."

                                                                         SLD   PSL    UP Catholic Other UD.KLD
(Intercept)/(Intercept)/(Intercept)/(Intercept)/(Intercept)/Suchoka.UD FALSE FALSE FALSE    FALSE FALSE   TRUE
New.Jobs/New.Jobs/New.Jobs/New.Jobs/New.Jobs/NA                        FALSE FALSE FALSE    FALSE FALSE  FALSE
Unemp/Unemp/Unemp/Unemp/Unemp/NA                                       FALSE FALSE FALSE    FALSE FALSE  FALSE
SOE/SOE/SOE/SOE/SOE/NA                                                 FALSE FALSE FALSE    FALSE FALSE  FALSE
Church.Att/Church.Att/Church.Att/Church.Att/Church.Att/NA              FALSE FALSE FALSE    FALSE FALSE  FALSE
Schooling/Schooling/Schooling/Schooling/Schooling/NA                   FALSE FALSE FALSE    FALSE FALSE  FALSE
Age/Age/Age/Age/Age/NA                                                 FALSE FALSE FALSE    FALSE FALSE  FALSE
Kwasn.SLD/Farmer/Zioqkowska.UP/NA/NA/NA                                 TRUE FALSE  TRUE    FALSE FALSE  FALSE
NA/Village/NA/NA/NA/NA                                                 FALSE FALSE FALSE    FALSE FALSE  FALSE
NA/Pawlak.PSL/NA/NA/NA/NA                                              FALSE  TRUE FALSE    FALSE FALSE  FALSE



multinomRob(): Grouped MNL Estimation
[1] "multinomMLE: hessian determinant: 8.32364479084216e+174"
[1] "multinomMLE: OPG determinant: -6.91388059961834e+232"
MNL LQD Fit: 45.34958 
MNL Estimates:
                                                                               SLD        PSL         UP   Catholic      Other
(Intercept)/(Intercept)/(Intercept)/(Intercept)/(Intercept)/Suchoka.UD  5.56470216  1.9283937  1.7559552 -6.8425136  2.8750465
New.Jobs/New.Jobs/New.Jobs/New.Jobs/New.Jobs/NA                        -6.38029005 -7.7064767 -3.7073963 -3.7837059 -5.2874834
Unemp/Unemp/Unemp/Unemp/Unemp/NA                                        0.56934823  0.8548404  0.7886248  1.6848328  0.6470373
SOE/SOE/SOE/SOE/SOE/NA                                                  0.38299624  0.4727334  0.7133912  0.3639858  0.8526141
Church.Att/Church.Att/Church.Att/Church.Att/Church.Att/NA              -2.04035701 -1.5831282 -1.8045749  2.3092508  0.1831089
Schooling/Schooling/Schooling/Schooling/Schooling/NA                   -0.38110595 -0.3512309 -0.2795588 -0.3455569 -0.4004014
Age/Age/Age/Age/Age/NA                                                  0.07937176  0.3992323  0.3623334  1.8365894  0.4991388
Kwasn.SLD/Farmer/Zioqkowska.UP/NA/NA/NA                                 0.52296642  1.9779548  0.9099635  0.0000000  0.0000000
NA/Village/NA/NA/NA/NA                                                  0.00000000  2.1906120  0.0000000  0.0000000  0.0000000
NA/Pawlak.PSL/NA/NA/NA/NA                                               0.00000000  0.9430359  0.0000000  0.0000000  0.0000000
                                                                         UD.KLD
(Intercept)/(Intercept)/(Intercept)/(Intercept)/(Intercept)/Suchoka.UD 0.554449
New.Jobs/New.Jobs/New.Jobs/New.Jobs/New.Jobs/NA                        0.000000
Unemp/Unemp/Unemp/Unemp/Unemp/NA                                       0.000000
SOE/SOE/SOE/SOE/SOE/NA                                                 0.000000
Church.Att/Church.Att/Church.Att/Church.Att/Church.Att/NA              0.000000
Schooling/Schooling/Schooling/Schooling/Schooling/NA                   0.000000
Age/Age/Age/Age/Age/NA                                                 0.000000
Kwasn.SLD/Farmer/Zioqkowska.UP/NA/NA/NA                                0.000000
NA/Village/NA/NA/NA/NA                                                 0.000000
NA/Pawlak.PSL/NA/NA/NA/NA                                              0.000000

MNL SEs:
                                                                              SLD        PSL         UP  Catholic      Other
(Intercept)/(Intercept)/(Intercept)/(Intercept)/(Intercept)/Suchoka.UD 2.14703642 3.65464929 2.32539814 3.8296067 2.02377999
New.Jobs/New.Jobs/New.Jobs/New.Jobs/New.Jobs/NA                        2.30884778 3.09051578 1.54883390 2.6504208 1.72559110
Unemp/Unemp/Unemp/Unemp/Unemp/NA                                       0.89485586 1.30009354 0.80165246 1.4986097 0.95648655
SOE/SOE/SOE/SOE/SOE/NA                                                 0.35653682 0.40899682 0.24045791 0.5652509 0.31320050
Church.Att/Church.Att/Church.Att/Church.Att/Church.Att/NA              0.41637138 0.53084736 0.41378032 0.6323097 0.40603024
Schooling/Schooling/Schooling/Schooling/Schooling/NA                   0.09580696 0.11460784 0.07528978 0.1517330 0.08688971
Age/Age/Age/Age/Age/NA                                                 0.36441460 0.64493494 0.37201770 0.6088072 0.32035252
Kwasn.SLD/Farmer/Zioqkowska.UP/NA/NA/NA                                0.04661507 0.92695473 0.07635412 0.0000000 0.00000000
NA/Village/NA/NA/NA/NA                                                 0.00000000 0.80455386 0.00000000 0.0000000 0.00000000
NA/Pawlak.PSL/NA/NA/NA/NA                                              0.00000000 0.07285661 0.00000000 0.0000000 0.00000000
                                                                           UD.KLD
(Intercept)/(Intercept)/(Intercept)/(Intercept)/(Intercept)/Suchoka.UD 0.08705206
New.Jobs/New.Jobs/New.Jobs/New.Jobs/New.Jobs/NA                        0.00000000
Unemp/Unemp/Unemp/Unemp/Unemp/NA                                       0.00000000
SOE/SOE/SOE/SOE/SOE/NA                                                 0.00000000
Church.Att/Church.Att/Church.Att/Church.Att/Church.Att/NA              0.00000000
Schooling/Schooling/Schooling/Schooling/Schooling/NA                   0.00000000
Age/Age/Age/Age/Age/NA                                                 0.00000000
Kwasn.SLD/Farmer/Zioqkowska.UP/NA/NA/NA                                0.00000000
NA/Village/NA/NA/NA/NA                                                 0.00000000
NA/Pawlak.PSL/NA/NA/NA/NA                                              0.00000000

multinomRob(): Using grouped MNL estimates as starting values.
multinomRob(): Starting Values 
                                                                               SLD        PSL         UP   Catholic      Other
(Intercept)/(Intercept)/(Intercept)/(Intercept)/(Intercept)/Suchoka.UD  5.56470216  1.9283937  1.7559552 -6.8425136  2.8750465
New.Jobs/New.Jobs/New.Jobs/New.Jobs/New.Jobs/NA                        -6.38029005 -7.7064767 -3.7073963 -3.7837059 -5.2874834
Unemp/Unemp/Unemp/Unemp/Unemp/NA                                        0.56934823  0.8548404  0.7886248  1.6848328  0.6470373
SOE/SOE/SOE/SOE/SOE/NA                                                  0.38299624  0.4727334  0.7133912  0.3639858  0.8526141
Church.Att/Church.Att/Church.Att/Church.Att/Church.Att/NA              -2.04035701 -1.5831282 -1.8045749  2.3092508  0.1831089
Schooling/Schooling/Schooling/Schooling/Schooling/NA                   -0.38110595 -0.3512309 -0.2795588 -0.3455569 -0.4004014
Age/Age/Age/Age/Age/NA                                                  0.07937176  0.3992323  0.3623334  1.8365894  0.4991388
Kwasn.SLD/Farmer/Zioqkowska.UP/NA/NA/NA                                 0.52296642  1.9779548  0.9099635  0.0000000  0.0000000
NA/Village/NA/NA/NA/NA                                                  0.00000000  2.1906120  0.0000000  0.0000000  0.0000000
NA/Pawlak.PSL/NA/NA/NA/NA                                               0.00000000  0.9430359  0.0000000  0.0000000  0.0000000
                                                                         UD.KLD
(Intercept)/(Intercept)/(Intercept)/(Intercept)/(Intercept)/Suchoka.UD 0.554449
New.Jobs/New.Jobs/New.Jobs/New.Jobs/New.Jobs/NA                        0.000000
Unemp/Unemp/Unemp/Unemp/Unemp/NA                                       0.000000
SOE/SOE/SOE/SOE/SOE/NA                                                 0.000000
Church.Att/Church.Att/Church.Att/Church.Att/Church.Att/NA              0.000000
Schooling/Schooling/Schooling/Schooling/Schooling/NA                   0.000000
Age/Age/Age/Age/Age/NA                                                 0.000000
Kwasn.SLD/Farmer/Zioqkowska.UP/NA/NA/NA                                0.000000
NA/Village/NA/NA/NA/NA                                                 0.000000
NA/Pawlak.PSL/NA/NA/NA/NA                                              0.000000
multinomRob(): starting fit = 45.34958 


Thu Mar  4 16:50:23 2004
Domains:
 -5.008232e+01   <=  X1   <=    6.121172e+01 
 -7.018319e+01   <=  X2   <=    5.742261e+01 
 -5.124134e+00   <=  X3   <=    6.262831e+00 
 -3.446966e+00   <=  X4   <=    4.212959e+00 
 -2.244393e+01   <=  X5   <=    1.836321e+01 
 -4.192165e+00   <=  X6   <=    3.429954e+00 
 -7.143459e-01   <=  X7   <=    8.730894e-01 
 -4.706698e+00   <=  X8   <=    5.752631e+00 
 -1.735554e+01   <=  X9   <=    2.121233e+01 
 -8.477124e+01   <=  X10  <=    6.935829e+01 
 -7.693564e+00   <=  X11  <=    9.403245e+00 
 -4.254600e+00   <=  X12  <=    5.200067e+00 
 -1.741441e+01   <=  X13  <=    1.424815e+01 
 -3.863540e+00   <=  X14  <=    3.161078e+00 
 -3.593091e+00   <=  X15  <=    4.391555e+00 
 -1.780159e+01   <=  X16  <=    2.175750e+01 
 -1.971551e+01   <=  X17  <=    2.409673e+01 
 -8.487323e+00   <=  X18  <=    1.037339e+01 
 -1.580360e+01   <=  X19  <=    1.931551e+01 
 -4.078136e+01   <=  X20  <=    3.336657e+01 
 -7.097623e+00   <=  X21  <=    8.674873e+00 
 -6.420521e+00   <=  X22  <=    7.847303e+00 
 -1.985032e+01   <=  X23  <=    1.624117e+01 
 -3.075147e+00   <=  X24  <=    2.516029e+00 
 -3.261000e+00   <=  X25  <=    3.985667e+00 
 -8.189671e+00   <=  X26  <=    1.000960e+01 
 -7.526765e+01   <=  X27  <=    6.158262e+01 
 -4.162077e+01   <=  X28  <=    3.405335e+01 
 -1.516350e+01   <=  X29  <=    1.853316e+01 
 -3.275872e+00   <=  X30  <=    4.003843e+00 
 -2.078326e+01   <=  X31  <=    2.540176e+01 
 -3.801125e+00   <=  X32  <=    3.110012e+00 
 -1.652930e+01   <=  X33  <=    2.020248e+01 
 -2.587542e+01   <=  X34  <=    3.162551e+01 
 -5.816232e+01   <=  X35  <=    4.758735e+01 
 -5.823336e+00   <=  X36  <=    7.117410e+00 
 -7.673527e+00   <=  X37  <=    9.378756e+00 
 -1.647980e+00   <=  X38  <=    2.014198e+00 
 -4.404416e+00   <=  X39  <=    3.603613e+00 
 -4.492249e+00   <=  X40  <=    5.490527e+00 
 -4.990041e+00   <=  X41  <=    6.098939e+00 

WARNING: Operator 6 (Multiple Point Simple Crossover) may only be started
WARNING: an even number of times.  I am increasing this operator by one.

WARNING: Operator 8 (Heuristic Crossover) may only be started
WARNING: an even number of times.  I am increasing this operator by one.

Data Type: Floating Point
Operators (code number, name, population) 
	(1) Cloning........................... 	122
	(2) Uniform Mutation.................. 	125
	(3) Boundary Mutation................. 	125
	(4) Non-Uniform Mutation.............. 	125
	(5) Polytope Crossover................ 	125
	(6) Multiple Point Simple Crossover... 	126
	(7) Whole Non-Uniform Mutation........ 	125
	(8) Heuristic Crossover............... 	126
	(9) Local-Minimum Crossover........... 	0

SOFT Maximum Number of Generations: 100
Maximum Nonchanging Generations: 10
Population size       : 1000
Convergence Tolerance: 1.000000e-07

Using the BFGS Derivative Based Optimizer on the Best Individual Each Generation.
Not Checking Gradients before Stoping
Using Out of Bounds Individuals.



Minimization Problem.



STARTING VALUES


The 2 best initial individuals are
 5.56	-6.38	 0.57	 0.38	-2.04	-0.38	 0.08	 0.52	 1.93	-7.71	 0.85	 0.47	-1.58	-0.35	 0.40	 1.98	 2.19	 0.94	 1.76	-3.71	 0.79	 0.71	-1.80	-0.28	 0.36	 0.91	-6.84	-3.78	 1.68	 0.36	 2.31	-0.35	 1.84	 2.88	-5.29	 0.65	 0.85	 0.18	-0.40	 0.50	 0.55	
fitness = 4.534958e+01

 7.63	46.17	-1.28	-2.32	15.81	-2.16	 0.63	 3.61	20.29	-22.12	 4.46	 3.19	 9.10	-3.72	-0.01	-3.72	20.18	 0.97	14.99	 0.62	 2.91	-0.16	 8.87	-1.48	-1.21	-6.55	14.54	27.25	-5.63	 0.17	-14.65	-2.75	-6.12	19.81	32.08	 0.35	-7.18	-0.64	-0.92	-1.17	 3.75	
fitness = 4.650732e+03


The worst fit of the population is: 1.227019e+40




Generation#	    Solution Value

      0 	4.534958e+01
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
      1 	4.495121e+01
XXX
XXX
      3 	4.495121e+01
XXX
XXX
XXX
XXX
XXX
XXX
      4 	4.483587e+01
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
      5 	4.471851e+01
XXX
XXX
XXX
      6 	4.455783e+01
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
      7 	4.454613e+01
XXX
XXX
XXX
      8 	4.453249e+01
      9 	4.445343e+01
XXX
XXX
XXX
     10 	4.435886e+01
XXX
     11 	4.430293e+01
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
     12 	4.416701e+01
     13 	4.416701e+01
     14 	4.415244e+01
     15 	4.414848e+01
     16 	4.414075e+01
XXX
XXX
     17 	4.412568e+01
XXX
     18 	4.412419e+01
XXX
XXX
     19 	4.411411e+01
XXX
XXX
     21 	4.411408e+01
XXX
     22 	4.411407e+01
XXX
XXX
XXX
XXX
     23 	4.410734e+01
XXX
XXX
     24 	4.409918e+01
XXX
     25 	4.408761e+01
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
     26 	4.406897e+01
XXX
     27 	4.406738e+01
XXX
     28 	4.406675e+01
XXX
     29 	4.406674e+01
XXX
     30 	4.406673e+01
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
     31 	4.406501e+01
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
     32 	4.405542e+01
XXX
XXX
     34 	4.405457e+01
XXX
     35 	4.405397e+01
XXX
     36 	4.405352e+01
XXX
     37 	4.404173e+01
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
     38 	4.403383e+01
XXX
     39 	4.402986e+01
XXX
     40 	4.402185e+01
XXX
     41 	4.402174e+01
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
     42 	4.401855e+01
XXX
     43 	4.401855e+01
XXX
     44 	4.401855e+01
XXX
     45 	4.401855e+01
XXX
XXX
XXX
XXX
XXX
XXX
     46 	4.401076e+01
XXX
     47 	4.401051e+01
XXX
     48 	4.401049e+01
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
     49 	4.400746e+01
XXX
     50 	4.400746e+01
XXX
XXX
     51 	4.397177e+01
XXX
XXX
XXX
     52 	4.386219e+01
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
     53 	4.379884e+01
     54 	4.375117e+01
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
     55 	4.356617e+01
XXX
     56 	4.356324e+01
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
     57 	4.352323e+01
     58 	4.351025e+01
     59 	4.350556e+01
     60 	4.349109e+01
     61 	4.348472e+01
     62 	4.348393e+01
     63 	4.347544e+01
XXX
XXX
XXX
     64 	4.346830e+01
XXX
XXX
     65 	4.345191e+01
     66 	4.344997e+01
     67 	4.344495e+01
     68 	4.343155e+01
XXX
XXX
     69 	4.342548e+01
     70 	4.342055e+01
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
     71 	4.340324e+01
     72 	4.340154e+01
     73 	4.339805e+01
     74 	4.339637e+01
     75 	4.339347e+01
     76 	4.339121e+01
     77 	4.338774e+01
     78 	4.338412e+01
XXX
XXX
     79 	4.338331e+01
     80 	4.338252e+01
     81 	4.337892e+01
     82 	4.337861e+01
     83 	4.337746e+01
     84 	4.337744e+01
XXX
XXX
XXX
XXX
XXX
XXX
     85 	4.337713e+01
     86 	4.337709e+01
     87 	4.337702e+01
     88 	4.337701e+01
     89 	4.337700e+01
     90 	4.337699e+01
     91 	4.337698e+01
     92 	4.337695e+01
     93 	4.337690e+01
     94 	4.337690e+01
     95 	4.337690e+01
     96 	4.337690e+01
     97 	4.337689e+01
XXX
XXX
XXX
XXX
XXX
XXX
XXX
     98 	4.337688e+01
    100 	4.337687e+01

Increasing Soft Maximum Generation Limit by 100 (MaxGenerations) to 200.
I'm doing this because the fitness is still improving.

    101 	4.337687e+01
    102 	4.337687e+01
    103 	4.337687e+01
    104 	4.337687e+01
    105 	4.337687e+01
    106 	4.337685e+01
XXX
    107 	4.337684e+01
    108 	4.337679e+01
XXX
XXX
    109 	4.337587e+01
XXX
XXX
XXX
XXX
    110 	4.337059e+01
XXX
XXX
    111 	4.337054e+01
    112 	4.337049e+01
    113 	4.337046e+01
    114 	4.337046e+01
    115 	4.336894e+01
    116 	4.336879e+01
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
    117 	4.336858e+01
XXX
XXX
XXX
    118 	4.336858e+01
XXX
XXX
XXX
XXX
    119 	4.336853e+01
    120 	4.336849e+01
    121 	4.336848e+01
    122 	4.336848e+01
XXX
XXX
XXX
XXX
XXX
    123 	4.336848e+01
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
    124 	4.336715e+01
XXX
    125 	4.336634e+01
    126 	4.336590e+01
XXX
    127 	4.336590e+01
    128 	4.336590e+01
    129 	4.336520e+01
    130 	4.336481e+01
    131 	4.336481e+01
    132 	4.336480e+01
    133 	4.336476e+01
    134 	4.336462e+01
    135 	4.336457e+01
    136 	4.336444e+01
    137 	4.336398e+01
    138 	4.336397e+01

WARNING: Turning Off MemoryMatrix because memory usage is too great.

    139 	4.336396e+01
    140 	4.336388e+01
    141 	4.336387e+01
    142 	4.336365e+01
    143 	4.336305e+01
    144 	4.336203e+01
    145 	4.336201e+01
    146 	4.336107e+01
    147 	4.336029e+01
    148 	4.335982e+01
    149 	4.335973e+01
    150 	4.335894e+01
    151 	4.335857e+01
    152 	4.335807e+01
    153 	4.335807e+01
    154 	4.335791e+01
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
    155 	4.335766e+01
    156 	4.335766e+01
XXX
XXX
XXX
XXX
    157 	4.335750e+01
    158 	4.335735e+01
    159 	4.335711e+01
XXX
XXX
XXX
XXX
XXX
XXX
XXX
    160 	4.335696e+01
XXX
XXX
XXX
    161 	4.335689e+01
XXX
XXX
XXX
    162 	4.335623e+01
    163 	4.335622e+01
    164 	4.335619e+01
    165 	4.335618e+01
    166 	4.335618e+01
    167 	4.335604e+01
    168 	4.335594e+01
    169 	4.335583e+01
    170 	4.335556e+01
    171 	4.335553e+01
    172 	4.335539e+01
    173 	4.335527e+01
    174 	4.335526e+01
    175 	4.335526e+01
    176 	4.335525e+01
    177 	4.335524e+01
    178 	4.335522e+01
    179 	4.335522e+01
    180 	4.335522e+01
    181 	4.335522e+01
    182 	4.335522e+01
    183 	4.335522e+01
    184 	4.335522e+01
    185 	4.335522e+01
    186 	4.335522e+01
    187 	4.335521e+01
    188 	4.335521e+01
XXX
XXX
XXX
XXX
XXX
    189 	4.335517e+01
    190 	4.335509e+01
    191 	4.335509e+01
    192 	4.335508e+01
    193 	4.335507e+01
    194 	4.335507e+01
    195 	4.335506e+01
    196 	4.335506e+01
    197 	4.335506e+01
    198 	4.335506e+01
    199 	4.335506e+01
    200 	4.335506e+01

Increasing Soft Maximum Generation Limit by 200 (MaxGenerations) to 400.
I'm doing this because the fitness is still improving.

    201 	4.335505e+01
    202 	4.335505e+01
    203 	4.335505e+01
    204 	4.335505e+01
    205 	4.335505e+01
    206 	4.335505e+01
    207 	4.335505e+01
    208 	4.335505e+01
XXX
XXX
XXX
XXX
    209 	4.335502e+01
    210 	4.335501e+01
    211 	4.335495e+01
    212 	4.335492e+01
    213 	4.335490e+01
    214 	4.335479e+01
    215 	4.335476e+01
    216 	4.335472e+01
    217 	4.335471e+01
    218 	4.335465e+01
    219 	4.335463e+01
    220 	4.335455e+01
    221 	4.335451e+01
    222 	4.335450e+01
    223 	4.335444e+01
    224 	4.335444e+01
    225 	4.335444e+01
    226 	4.335444e+01
    227 	4.335444e+01
    228 	4.335444e+01
    229 	4.335444e+01
    230 	4.335444e+01
    231 	4.335444e+01
    232 	4.335444e+01
    233 	4.335444e+01
    234 	4.335444e+01
    235 	4.335444e+01
    236 	4.335444e+01
    237 	4.335444e+01
    238 	4.335444e+01
    239 	4.335443e+01
    240 	4.335443e+01
    241 	4.335442e+01
    242 	4.335442e+01
    243 	4.335442e+01
    244 	4.335442e+01
    245 	4.335442e+01
    246 	4.335442e+01
    247 	4.335442e+01
    248 	4.335442e+01
    249 	4.335442e+01
    250 	4.335442e+01
    251 	4.335442e+01
    252 	4.335442e+01
    253 	4.335442e+01

Soft Generation Wait Limit Hit.
No Improvement in 10 Generations

Best Fit Found at Generation 253
Fit Value = 4.335442e+01


Parameters at the Solution (value, gradient):

 X[ 1] :	5.754089e+00	G[ 1] :	-1.114518e+02
 X[ 2] :	-6.378191e+00	G[ 2] :	-3.842926e+00
 X[ 3] :	5.693368e-01	G[ 3] :	-1.346809e+01
 X[ 4] :	3.831981e-01	G[ 4] :	-2.937873e+01
 X[ 5] :	-2.351269e+00	G[ 5] :	-8.236780e+01
 X[ 6] :	-3.839953e-01	G[ 6] :	-1.284969e+03
 X[ 7] :	7.819353e-02	G[ 7] :	-4.834368e+02
 X[ 8] :	5.658545e-01	G[ 8] :	0.000000e+00
 X[ 9] :	1.927910e+00	G[ 9] :	3.658597e+01
 X[10] :	-7.706419e+00	G[10] :	3.739209e-01
 X[11] :	8.533632e-01	G[11] :	9.434852e-01
 X[12] :	4.512534e-01	G[12] :	4.997135e+00
 X[13] :	-1.583214e+00	G[13] :	4.927002e+00
 X[14] :	-3.489794e-01	G[14] :	-9.546280e+01
 X[15] :	4.060728e-01	G[15] :	4.189190e+01
 X[16] :	1.884470e+00	G[16] :	2.977418e+01
 X[17] :	2.114824e+00	G[17] :	-2.183413e+01
 X[18] :	9.537247e-01	G[18] :	0.000000e+00
 X[19] :	1.756289e+00	G[19] :	2.612690e+02
 X[20] :	-3.676550e+00	G[20] :	4.500009e+01
 X[21] :	6.924822e-01	G[21] :	1.195641e+01
 X[22] :	7.111356e-01	G[22] :	2.396374e+02
 X[23] :	-1.804689e+00	G[23] :	1.177286e+02
 X[24] :	-2.813908e-01	G[24] :	3.325639e+03
 X[25] :	3.618087e-01	G[25] :	1.195478e+03
 X[26] :	7.750130e-01	G[26] :	0.000000e+00
 X[27] :	-6.842894e+00	G[27] :	3.384005e+00
 X[28] :	-3.783658e+00	G[28] :	-1.578616e+00
 X[29] :	1.685790e+00	G[29] :	1.595918e+00
 X[30] :	-3.350631e-01	G[30] :	-8.119422e+00
 X[31] :	2.308825e+00	G[31] :	4.933464e+00
 X[32] :	-3.458717e-01	G[32] :	1.018606e+01
 X[33] :	1.836247e+00	G[33] :	1.234167e+02
 X[34] :	2.875250e+00	G[34] :	-7.511451e+01
 X[35] :	-5.287481e+00	G[35] :	-2.341630e+01
 X[36] :	6.477857e-01	G[36] :	3.002903e+00
 X[37] :	8.563832e-01	G[37] :	-1.213630e+02
 X[38] :	1.768028e-01	G[38] :	-1.494446e+01
 X[39] :	-4.004084e-01	G[39] :	3.285060e+01
 X[40] :	5.006597e-01	G[40] :	5.477607e+01
 X[41] :	-1.692858e-01	G[41] :	0.000000e+00

Solution Found Generation 253
Number of Generations Run 253


Thu Mar  4 17:02:01 2004
Total run time : 0 hours 11 minutes and 38 seconds


LQD Results:
                                                                               SLD        PSL         UP   Catholic      Other
(Intercept)/(Intercept)/(Intercept)/(Intercept)/(Intercept)/Suchoka.UD  5.75408895  1.9279097  1.7562887 -6.8428942  2.8752501
New.Jobs/New.Jobs/New.Jobs/New.Jobs/New.Jobs/NA                        -6.37819132 -7.7064185 -3.6765498 -3.7836582 -5.2874815
Unemp/Unemp/Unemp/Unemp/Unemp/NA                                        0.56933678  0.8533632  0.6924822  1.6857904  0.6477857
SOE/SOE/SOE/SOE/SOE/NA                                                  0.38319811  0.4512534  0.7111356 -0.3350631  0.8563832
Church.Att/Church.Att/Church.Att/Church.Att/Church.Att/NA              -2.35126935 -1.5832140 -1.8046893  2.3088249  0.1768028
Schooling/Schooling/Schooling/Schooling/Schooling/NA                   -0.38399528 -0.3489794 -0.2813908 -0.3458717 -0.4004084
Age/Age/Age/Age/Age/NA                                                  0.07819353  0.4060728  0.3618087  1.8362465  0.5006597
Kwasn.SLD/Farmer/Zioqkowska.UP/NA/NA/NA                                 0.56585448  1.8844699  0.7750130  0.0000000  0.0000000
NA/Village/NA/NA/NA/NA                                                  0.00000000  2.1148236  0.0000000  0.0000000  0.0000000
NA/Pawlak.PSL/NA/NA/NA/NA                                               0.00000000  0.9537247  0.0000000  0.0000000  0.0000000
                                                                           UD.KLD
(Intercept)/(Intercept)/(Intercept)/(Intercept)/(Intercept)/Suchoka.UD -0.1692858
New.Jobs/New.Jobs/New.Jobs/New.Jobs/New.Jobs/NA                         0.0000000
Unemp/Unemp/Unemp/Unemp/Unemp/NA                                        0.0000000
SOE/SOE/SOE/SOE/SOE/NA                                                  0.0000000
Church.Att/Church.Att/Church.Att/Church.Att/Church.Att/NA               0.0000000
Schooling/Schooling/Schooling/Schooling/Schooling/NA                    0.0000000
Age/Age/Age/Age/Age/NA                                                  0.0000000
Kwasn.SLD/Farmer/Zioqkowska.UP/NA/NA/NA                                 0.0000000
NA/Village/NA/NA/NA/NA                                                  0.0000000
NA/Pawlak.PSL/NA/NA/NA/NA                                               0.0000000

LQD sigma: 43.35442 


(multinomTanh):
Tanh Estimates
                                                                               SLD        PSL         UP   Catholic      Other
(Intercept)/(Intercept)/(Intercept)/(Intercept)/(Intercept)/Suchoka.UD  5.76281359  1.9318264  1.7575815 -6.8498618  2.8700406
New.Jobs/New.Jobs/New.Jobs/New.Jobs/New.Jobs/NA                        -6.38161390 -7.6904903 -3.6701397 -3.5706507 -5.2926677
Unemp/Unemp/Unemp/Unemp/Unemp/NA                                        0.55755660  0.8602874  0.7163444  1.8011899  0.6483360
SOE/SOE/SOE/SOE/SOE/NA                                                  0.39045127  0.4559892  0.7150949 -0.1555647  0.8559880
Church.Att/Church.Att/Church.Att/Church.Att/Church.Att/NA              -2.30890827 -1.5775708 -1.8008069  2.3497920  0.1780828
Schooling/Schooling/Schooling/Schooling/Schooling/NA                   -0.38475277 -0.3498159 -0.2807362 -0.3472559 -0.4000282
Age/Age/Age/Age/Age/NA                                                  0.07312998  0.4043774  0.3593593  1.8245505  0.5006840
Kwasn.SLD/Farmer/Zioqkowska.UP/NA/NA/NA                                 0.52416632  1.8996130  0.9095500  0.0000000  0.0000000
NA/Village/NA/NA/NA/NA                                                  0.00000000  2.1202716  0.0000000  0.0000000  0.0000000
NA/Pawlak.PSL/NA/NA/NA/NA                                               0.00000000  0.9433627  0.0000000  0.0000000  0.0000000
                                                                          UD.KLD
(Intercept)/(Intercept)/(Intercept)/(Intercept)/(Intercept)/Suchoka.UD 0.5532653
New.Jobs/New.Jobs/New.Jobs/New.Jobs/New.Jobs/NA                        0.0000000
Unemp/Unemp/Unemp/Unemp/Unemp/NA                                       0.0000000
SOE/SOE/SOE/SOE/SOE/NA                                                 0.0000000
Church.Att/Church.Att/Church.Att/Church.Att/Church.Att/NA              0.0000000
Schooling/Schooling/Schooling/Schooling/Schooling/NA                   0.0000000
Age/Age/Age/Age/Age/NA                                                 0.0000000
Kwasn.SLD/Farmer/Zioqkowska.UP/NA/NA/NA                                0.0000000
NA/Village/NA/NA/NA/NA                                                 0.0000000
NA/Pawlak.PSL/NA/NA/NA/NA                                              0.0000000

Tanh Sandwich SEs
                                                                              SLD        PSL         UP  Catholic      Other
(Intercept)/(Intercept)/(Intercept)/(Intercept)/(Intercept)/Suchoka.UD 2.08510493 3.67229020 2.33721238 5.4468754 2.01984981
New.Jobs/New.Jobs/New.Jobs/New.Jobs/New.Jobs/NA                        2.24472920 3.10346268 1.38618638 3.8860846 1.69962969
Unemp/Unemp/Unemp/Unemp/Unemp/NA                                       0.83939329 1.27773893 0.80134752 1.8491326 0.95112434
SOE/SOE/SOE/SOE/SOE/NA                                                 0.34705329 0.41134956 0.24562522 0.9571347 0.31244175
Church.Att/Church.Att/Church.Att/Church.Att/Church.Att/NA              0.40177424 0.51887707 0.42467403 0.7938496 0.39985160
Schooling/Schooling/Schooling/Schooling/Schooling/NA                   0.09120428 0.11382597 0.07525722 0.1629559 0.08522077
Age/Age/Age/Age/Age/NA                                                 0.35719290 0.64778443 0.35696591 0.9083497 0.32011721
Kwasn.SLD/Farmer/Zioqkowska.UP/NA/NA/NA                                0.04718871 0.92033691 0.07985506 0.0000000 0.00000000
NA/Village/NA/NA/NA/NA                                                 0.00000000 0.79291208 0.00000000 0.0000000 0.00000000
NA/Pawlak.PSL/NA/NA/NA/NA                                              0.00000000 0.07270573 0.00000000 0.0000000 0.00000000
                                                                         UD.KLD
(Intercept)/(Intercept)/(Intercept)/(Intercept)/(Intercept)/Suchoka.UD 0.085175
New.Jobs/New.Jobs/New.Jobs/New.Jobs/New.Jobs/NA                        0.000000
Unemp/Unemp/Unemp/Unemp/Unemp/NA                                       0.000000
SOE/SOE/SOE/SOE/SOE/NA                                                 0.000000
Church.Att/Church.Att/Church.Att/Church.Att/Church.Att/NA              0.000000
Schooling/Schooling/Schooling/Schooling/Schooling/NA                   0.000000
Age/Age/Age/Age/Age/NA                                                 0.000000
Kwasn.SLD/Farmer/Zioqkowska.UP/NA/NA/NA                                0.000000
NA/Village/NA/NA/NA/NA                                                 0.000000
NA/Pawlak.PSL/NA/NA/NA/NA                                              0.000000

TANH sigma: 37.25687 


**********************************

Iteration: 1 
using best non-genoud LQD sigma. Provided by Tanh.


Thu Mar  4 17:02:03 2004
Domains:
 -2.305125e+01   <=  X1   <=    3.457688e+01 
 -3.828968e+01   <=  X2   <=    2.552646e+01 
 -2.230226e+00   <=  X3   <=    3.345340e+00 
 -1.561805e+00   <=  X4   <=    2.342708e+00 
 -1.385345e+01   <=  X5   <=    9.235633e+00 
 -2.308517e+00   <=  X6   <=    1.539011e+00 
 -2.925199e-01   <=  X7   <=    4.387799e-01 
 -2.096665e+00   <=  X8   <=    3.144998e+00 
 -7.727305e+00   <=  X9   <=    1.159096e+01 
 -4.614294e+01   <=  X10  <=    3.076196e+01 
 -3.441150e+00   <=  X11  <=    5.161724e+00 
 -1.823957e+00   <=  X12  <=    2.735935e+00 
 -9.465425e+00   <=  X13  <=    6.310283e+00 
 -2.098896e+00   <=  X14  <=    1.399264e+00 
 -1.617510e+00   <=  X15  <=    2.426265e+00 
 -7.598452e+00   <=  X16  <=    1.139768e+01 
 -8.481086e+00   <=  X17  <=    1.272163e+01 
 -3.773451e+00   <=  X18  <=    5.660176e+00 
 -7.030326e+00   <=  X19  <=    1.054549e+01 
 -2.202084e+01   <=  X20  <=    1.468056e+01 
 -2.865378e+00   <=  X21  <=    4.298067e+00 
 -2.860380e+00   <=  X22  <=    4.290570e+00 
 -1.080484e+01   <=  X23  <=    7.203228e+00 
 -1.684417e+00   <=  X24  <=    1.122945e+00 
 -1.437437e+00   <=  X25  <=    2.156156e+00 
 -3.638200e+00   <=  X26  <=    5.457300e+00 
 -4.109917e+01   <=  X27  <=    2.739945e+01 
 -2.142390e+01   <=  X28  <=    1.428260e+01 
 -7.204760e+00   <=  X29  <=    1.080714e+01 
 -9.333883e-01   <=  X30  <=    6.222589e-01 
 -9.399168e+00   <=  X31  <=    1.409875e+01 
 -2.083536e+00   <=  X32  <=    1.389024e+00 
 -7.298202e+00   <=  X33  <=    1.094730e+01 
 -1.148016e+01   <=  X34  <=    1.722024e+01 
 -3.175601e+01   <=  X35  <=    2.117067e+01 
 -2.593344e+00   <=  X36  <=    3.890016e+00 
 -3.423952e+00   <=  X37  <=    5.135928e+00 
 -7.123311e-01   <=  X38  <=    1.068497e+00 
 -2.400169e+00   <=  X39  <=    1.600113e+00 
 -2.002736e+00   <=  X40  <=    3.004104e+00 
 -2.213061e+00   <=  X41  <=    3.319592e+00 

WARNING: Operator 6 (Multiple Point Simple Crossover) may only be started
WARNING: an even number of times.  I am increasing this operator by one.

WARNING: Operator 8 (Heuristic Crossover) may only be started
WARNING: an even number of times.  I am increasing this operator by one.

Data Type: Floating Point
Operators (code number, name, population) 
	(1) Cloning........................... 	122
	(2) Uniform Mutation.................. 	125
	(3) Boundary Mutation................. 	125
	(4) Non-Uniform Mutation.............. 	125
	(5) Polytope Crossover................ 	125
	(6) Multiple Point Simple Crossover... 	126
	(7) Whole Non-Uniform Mutation........ 	125
	(8) Heuristic Crossover............... 	126
	(9) Local-Minimum Crossover........... 	0

SOFT Maximum Number of Generations: 100
Maximum Nonchanging Generations: 10
Population size       : 1000
Convergence Tolerance: 1.000000e-07

Using the BFGS Derivative Based Optimizer on the Best Individual Each Generation.
Not Checking Gradients before Stoping
Using Out of Bounds Individuals.



Minimization Problem.



STARTING VALUES


The 2 best initial individuals are
 5.76	-6.38	 0.56	 0.39	-2.31	-0.38	 0.07	 0.52	 1.93	-7.69	 0.86	 0.46	-1.58	-0.35	 0.40	 1.90	 2.12	 0.94	 1.76	-3.67	 0.72	 0.72	-1.80	-0.28	 0.36	 0.91	-6.85	-3.57	 1.80	-0.16	 2.35	-0.35	 1.82	 2.87	-5.29	 0.65	 0.86	 0.18	-0.40	 0.50	 0.55	
fitness = 4.482327e+01

17.23	-16.72	-0.08	 0.27	-7.77	-1.93	 0.43	 1.31	 9.68	 2.29	-2.88	-1.76	-3.37	-1.26	 0.86	-0.91	-2.66	-2.49	-4.85	-15.91	 1.89	 1.60	 3.28	-0.43	 0.78	 4.39	 6.62	-18.28	 5.27	-0.07	13.98	 0.14	-5.20	-7.00	-14.97	 0.53	 0.61	-0.43	-0.29	 1.37	-0.99	
fitness = 1.062474e+03


The worst fit of the population is: 4.048143e+19




Generation#	    Solution Value

      0 	4.482327e+01
XXX
XXX
XXX
XXX
XXX
      1 	4.467580e+01
XXX
XXX
XXX
      2 	4.416916e+01
      3 	4.411264e+01
XXX
XXX
      4 	4.377667e+01
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
      5 	4.369399e+01
XXX
XXX
XXX
XXX
XXX
XXX
      6 	4.367766e+01
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
      7 	4.365265e+01
      8 	4.359245e+01
      9 	4.357677e+01
XXX
XXX
XXX
XXX
XXX
     10 	4.354089e+01
XXX
     11 	4.353575e+01
XXX
XXX
XXX
XXX
     12 	4.352633e+01
XXX
     13 	4.352295e+01
XXX
     14 	4.352210e+01
XXX
     15 	4.352200e+01
XXX
     16 	4.351751e+01
XXX
     17 	4.351303e+01
XXX
     18 	4.351240e+01
XXX
     19 	4.351231e+01
XXX
     20 	4.351231e+01
XXX
     21 	4.351231e+01
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
     22 	4.350803e+01
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
     23 	4.349326e+01
     24 	4.347737e+01
     25 	4.347668e+01
XXX
     26 	4.347202e+01
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
     27 	4.347037e+01
     28 	4.346791e+01
     29 	4.346616e+01
     30 	4.346402e+01
     31 	4.346200e+01
     32 	4.346111e+01
     33 	4.346072e+01
     34 	4.345865e+01
     35 	4.345496e+01
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
     37 	4.345326e+01
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
     38 	4.344507e+01
XXX
     39 	4.343704e+01
XXX
     40 	4.343500e+01
XXX
     41 	4.343304e+01
XXX
XXX
     42 	4.342883e+01
XXX
     43 	4.342864e+01
XXX
     44 	4.342822e+01
XXX
     45 	4.342797e+01
XXX
     46 	4.342794e+01
XXX
     47 	4.342791e+01
XXX
     48 	4.342791e+01
XXX
     49 	4.342787e+01
XXX
     50 	4.342787e+01
XXX
     51 	4.342619e+01
XXX
XXX
     52 	4.342602e+01
XXX
XXX
     53 	4.342595e+01
XXX
XXX
     54 	4.342594e+01
XXX
XXX
XXX
XXX
XXX
XXX
XXX
     55 	4.342273e+01
XXX
XXX
XXX
XXX
     57 	4.341824e+01
XXX
XXX
XXX
     58 	4.341432e+01
     59 	4.341250e+01
XXX
XXX
     60 	4.340913e+01
     61 	4.340851e+01
     62 	4.340531e+01
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
     63 	4.340456e+01
XXX
XXX
XXX
XXX
XXX
     64 	4.340311e+01
     65 	4.340293e+01
     66 	4.340292e+01
     67 	4.340292e+01
     68 	4.340223e+01
     69 	4.340187e+01
     70 	4.339977e+01
     71 	4.339860e+01
XXX
XXX
     72 	4.339661e+01
     73 	4.339457e+01
     74 	4.339343e+01
     75 	4.339198e+01
XXX
XXX
XXX
XXX
     76 	4.338784e+01
     77 	4.338488e+01
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
     78 	4.338379e+01
     79 	4.338145e+01
     80 	4.338108e+01
     81 	4.338084e+01
     82 	4.338061e+01
     83 	4.338057e+01
     84 	4.338034e+01
     85 	4.338017e+01
     86 	4.338005e+01
     87 	4.337994e+01
     88 	4.337993e+01
     89 	4.337990e+01
     90 	4.337988e+01
     91 	4.337986e+01
XXX
     92 	4.337986e+01
XXX
     93 	4.337985e+01
XXX
     94 	4.337985e+01
XXX
     95 	4.337985e+01
XXX
     96 	4.337985e+01
XXX
     97 	4.337984e+01
XXX
     98 	4.337984e+01
XXX
     99 	4.337984e+01
XXX
    100 	4.337984e+01

Increasing Soft Maximum Generation Limit by 100 (MaxGenerations) to 200.
I'm doing this because the fitness is still improving.

XXX
    101 	4.337984e+01
XXX
    102 	4.337984e+01
XXX
    103 	4.337984e+01
XXX
    104 	4.337984e+01
XXX
    105 	4.337984e+01
XXX
    106 	4.337984e+01
XXX
    107 	4.337984e+01
XXX
    108 	4.337984e+01
XXX
    109 	4.337984e+01
XXX
    110 	4.337984e+01
XXX
    111 	4.337984e+01
XXX
    112 	4.337984e+01
XXX
    113 	4.337984e+01
XXX
    114 	4.337984e+01
XXX
    115 	4.337983e+01
    116 	4.337983e+01
XXX
    117 	4.337983e+01
XXX
    118 	4.337983e+01
XXX
    119 	4.337983e+01
XXX
    120 	4.337983e+01
XXX
    121 	4.337983e+01
XXX
    122 	4.337983e+01
XXX
    123 	4.337983e+01
XXX
    124 	4.337983e+01
XXX
    125 	4.337983e+01
XXX
    126 	4.337983e+01
XXX
    127 	4.337983e+01
XXX
    128 	4.337983e+01
XXX
    129 	4.337983e+01
XXX
    130 	4.337983e+01

Soft Generation Wait Limit Hit.
No Improvement in 10 Generations

Best Fit Found at Generation 130
Fit Value = 4.337983e+01


Parameters at the Solution (value, gradient):

 X[ 1] :	5.762860e+00	G[ 1] :	6.301868e+00
 X[ 2] :	-6.382835e+00	G[ 2] :	4.148315e-01
 X[ 3] :	5.297067e-01	G[ 3] :	9.785461e-01
 X[ 4] :	3.899018e-01	G[ 4] :	1.276023e+00
 X[ 5] :	-2.308863e+00	G[ 5] :	-6.358768e+01
 X[ 6] :	-3.847867e-01	G[ 6] :	-3.699970e+03
 X[ 7] :	7.333959e-02	G[ 7] :	-1.082586e+03
 X[ 8] :	6.003176e-01	G[ 8] :	0.000000e+00
 X[ 9] :	1.931603e+00	G[ 9] :	-1.042537e+02
 X[10] :	-7.690766e+00	G[10] :	1.322744e-01
 X[11] :	8.590907e-01	G[11] :	1.442774e-01
 X[12] :	4.558730e-01	G[12] :	-2.672967e+01
 X[13] :	-1.577675e+00	G[13] :	-1.449619e+02
 X[14] :	-3.496203e-01	G[14] :	-4.439981e+02
 X[15] :	4.043564e-01	G[15] :	-7.881959e+02
 X[16] :	1.902135e+00	G[16] :	-6.415666e-01
 X[17] :	2.120262e+00	G[17] :	-1.208219e+00
 X[18] :	9.557156e-01	G[18] :	0.000000e+00
 X[19] :	1.757725e+00	G[19] :	-7.289919e+01
 X[20] :	-3.670112e+00	G[20] :	3.091734e+00
 X[21] :	7.323209e-01	G[21] :	7.057264e+00
 X[22] :	7.155718e-01	G[22] :	-4.626417e+01
 X[23] :	-1.800748e+00	G[23] :	-1.485978e+01
 X[24] :	-2.805388e-01	G[24] :	-2.617016e+02
 X[25] :	3.594729e-01	G[25] :	-5.146988e+02
 X[26] :	7.849726e-01	G[26] :	0.000000e+00
 X[27] :	-6.849928e+00	G[27] :	8.719930e+00
 X[28] :	-3.570841e+00	G[28] :	-6.326818e-01
 X[29] :	1.801188e+00	G[29] :	8.541307e+00
 X[30] :	-1.593964e-01	G[30] :	7.455591e-01
 X[31] :	2.349787e+00	G[31] :	-1.389652e+01
 X[32] :	-3.474476e-01	G[32] :	1.195713e+02
 X[33] :	1.824439e+00	G[33] :	3.097177e+01
 X[34] :	2.870047e+00	G[34] :	4.476955e+00
 X[35] :	-5.292734e+00	G[35] :	1.173849e+00
 X[36] :	6.951229e-01	G[36] :	-2.696939e+01
 X[37] :	8.597068e-01	G[37] :	-8.118753e+01
 X[38] :	1.780703e-01	G[38] :	3.559766e+01
 X[39] :	-4.003952e-01	G[39] :	-2.420762e+00
 X[40] :	5.004511e-01	G[40] :	2.341596e+01
 X[41] :	-1.274926e-01	G[41] :	0.000000e+00

Solution Found Generation 130
Number of Generations Run 130


Thu Mar  4 17:07:25 2004
Total run time : 0 hours 5 minutes and 22 seconds

(iter 1 ): LQD Results:
                                                                              SLD        PSL         UP   Catholic      Other
(Intercept)/(Intercept)/(Intercept)/(Intercept)/(Intercept)/Suchoka.UD  5.7628595  1.9316030  1.7577247 -6.8499278  2.8700472
New.Jobs/New.Jobs/New.Jobs/New.Jobs/New.Jobs/NA                        -6.3828350 -7.6907660 -3.6701118 -3.5708410 -5.2927337
Unemp/Unemp/Unemp/Unemp/Unemp/NA                                        0.5297067  0.8590907  0.7323209  1.8011885  0.6951229
SOE/SOE/SOE/SOE/SOE/NA                                                  0.3899018  0.4558730  0.7155718 -0.1593964  0.8597068
Church.Att/Church.Att/Church.Att/Church.Att/Church.Att/NA              -2.3088629 -1.5776755 -1.8007484  2.3497872  0.1780703
Schooling/Schooling/Schooling/Schooling/Schooling/NA                   -0.3847867 -0.3496203 -0.2805388 -0.3474476 -0.4003952
Age/Age/Age/Age/Age/NA                                                  0.0733396  0.4043564  0.3594729  1.8244385  0.5004511
Kwasn.SLD/Farmer/Zioqkowska.UP/NA/NA/NA                                 0.6003176  1.9021352  0.7849726  0.0000000  0.0000000
NA/Village/NA/NA/NA/NA                                                  0.0000000  2.1202619  0.0000000  0.0000000  0.0000000
NA/Pawlak.PSL/NA/NA/NA/NA                                               0.0000000  0.9557156  0.0000000  0.0000000  0.0000000
                                                                           UD.KLD
(Intercept)/(Intercept)/(Intercept)/(Intercept)/(Intercept)/Suchoka.UD -0.1274926
New.Jobs/New.Jobs/New.Jobs/New.Jobs/New.Jobs/NA                         0.0000000
Unemp/Unemp/Unemp/Unemp/Unemp/NA                                        0.0000000
SOE/SOE/SOE/SOE/SOE/NA                                                  0.0000000
Church.Att/Church.Att/Church.Att/Church.Att/Church.Att/NA               0.0000000
Schooling/Schooling/Schooling/Schooling/Schooling/NA                    0.0000000
Age/Age/Age/Age/Age/NA                                                  0.0000000
Kwasn.SLD/Farmer/Zioqkowska.UP/NA/NA/NA                                 0.0000000
NA/Village/NA/NA/NA/NA                                                  0.0000000
NA/Pawlak.PSL/NA/NA/NA/NA                                               0.0000000

(iter 1 ): LQD sigma: 43.37983 

WARNING: New fit is worse.  Stopping.
It may be a good idea to restart with a larger GENOUD population.
> 
> summary(muljj, weights=TRUE);

Choice 1 : SLD Estimates and SE:
                  Est SE.Sand  SE.OPG SE.Hess t.val.Sand
(Intercept)    5.7600  2.0900 0.13700 0.07190      2.760
New.Jobs      -6.3800  2.2400 0.12200 0.06360     -2.840
Unemp          0.5580  0.8390 0.06870 0.02870      0.664
SOE            0.3900  0.3470 0.05640 0.01000      1.130
Church.Att    -2.3100  0.4020 0.02920 0.01320     -5.750
Schooling     -0.3850  0.0912 0.00587 0.00281     -4.220
Age            0.0731  0.3570 0.02810 0.01240      0.205
Kwasn.SLD      0.5240  0.0472 0.13200 0.00567     11.100
NA             0.0000  0.0000 0.00000 0.00000        NaN
NA             0.0000  0.0000 0.00000 0.00000        NaN


Choice 2 : PSL Estimates and SE:
                 Est SE.Sand  SE.OPG SE.Hess t.val.Sand
(Intercept)    1.930  3.6700 0.10200 0.08990      0.526
New.Jobs      -7.690  3.1000 0.11700 0.07660     -2.480
Unemp          0.860  1.2800 0.06640 0.03430      0.673
SOE            0.456  0.4110 0.06190 0.01230      1.110
Church.Att    -1.580  0.5190 0.02670 0.01600     -3.040
Schooling     -0.350  0.1140 0.00684 0.00315     -3.070
Age            0.404  0.6480 0.02760 0.01710      0.624
Farmer         1.900  0.9200 0.04240 0.02130      2.060
Village        2.120  0.7930 0.02150 0.01620      2.670
Pawlak.PSL     0.943  0.0727 0.11800 0.00500     13.000


Choice 3 : UP Estimates and SE:
                 Est SE.Sand SE.OPG SE.Hess t.val.Sand
(Intercept)    1.760  2.3400 0.1670 0.09650      0.752
New.Jobs      -3.670  1.3900 0.1500 0.08560     -2.650
Unemp          0.716  0.8010 0.1180 0.03910      0.894
SOE            0.715  0.2460 0.1210 0.01350      2.910
Church.Att    -1.800  0.4250 0.0578 0.01790     -4.240
Schooling     -0.281  0.0753 0.0110 0.00375     -3.730
Age            0.359  0.3570 0.0366 0.01650      1.010
Zioqkowska.UP  0.910  0.0799 0.0192 0.00509     11.400
NA             0.000  0.0000 0.0000 0.00000        NaN
NA             0.000  0.0000 0.0000 0.00000        NaN


Choice 4 : Catholic Estimates and SE:
                 Est SE.Sand  SE.OPG SE.Hess t.val.Sand
(Intercept)   -6.850   5.450 0.14800 0.11100     -1.260
New.Jobs      -3.570   3.890 0.13100 0.09320     -0.919
Unemp          1.800   1.850 0.09320 0.04230      0.974
SOE           -0.156   0.957 0.02960 0.01630     -0.163
Church.Att     2.350   0.794 0.04040 0.01720      2.960
Schooling     -0.347   0.163 0.00659 0.00396     -2.130
Age            1.820   0.908 0.02590 0.01940      2.010
NA             0.000   0.000 0.00000 0.00000        NaN
NA             0.000   0.000 0.00000 0.00000        NaN
NA             0.000   0.000 0.00000 0.00000        NaN


Choice 5 : Other Estimates and SE:
                 Est SE.Sand  SE.OPG SE.Hess t.val.Sand
(Intercept)    2.870  2.0200 0.08340 0.06550      1.420
New.Jobs      -5.290  1.7000 0.12400 0.05630     -3.110
Unemp          0.648  0.9510 0.06110 0.02630      0.682
SOE            0.856  0.3120 0.06460 0.00904      2.740
Church.Att     0.178  0.4000 0.02430 0.01150      0.445
Schooling     -0.400  0.0852 0.00568 0.00260     -4.690
Age            0.501  0.3200 0.01920 0.01130      1.560
NA             0.000  0.0000 0.00000 0.00000        NaN
NA             0.000  0.0000 0.00000 0.00000        NaN
NA             0.000  0.0000 0.00000 0.00000        NaN


Choice 6 : UD.KLD Estimates and SE:
                Est SE.Sand SE.OPG SE.Hess t.val.Sand
Suchoka.UD    0.553  0.0852 0.0427 0.00385        6.5
NA            0.000  0.0000 0.0000 0.00000        NaN
NA            0.000  0.0000 0.0000 0.00000        NaN
NA            0.000  0.0000 0.0000 0.00000        NaN
NA            0.000  0.0000 0.0000 0.00000        NaN
NA            0.000  0.0000 0.0000 0.00000        NaN
NA            0.000  0.0000 0.0000 0.00000        NaN
NA            0.000  0.0000 0.0000 0.00000        NaN
NA            0.000  0.0000 0.0000 0.00000        NaN
NA            0.000  0.0000 0.0000 0.00000        NaN


LQD sigma: 43.35442 
TANH sigma: 37.25687 

Number of Observations: 49 
Number of observations with at least one zero weight: 0 
Number of zero weights: 0 
TANH: weights
   name weights:SLD weights:PSL weights:UP weights:Catholic weights:Other
1     1     1.00000    1.000000          1         0.589126             1
2     2     1.00000    1.000000          1         1.000000             1
3     3     1.00000    0.871044          1         0.281952             1
4     4     1.00000    1.000000          1         1.000000             1
5     5     0.56820    1.000000          1         1.000000             1
6     6     1.00000    1.000000          1         1.000000             1
7     7     1.00000    1.000000          1         1.000000             1
8     8     1.00000    1.000000          1         1.000000             1
9     9     1.00000    1.000000          1         1.000000             1
10   10     1.00000    1.000000          1         1.000000             1
11   11     1.00000    1.000000          1         1.000000             1
12   12     1.00000    1.000000          1         1.000000             1
13   13     0.71112    1.000000          1         1.000000             1
14   14     1.00000    1.000000          1         1.000000             1
15   15     1.00000    0.943388          1         1.000000             1
16   16     1.00000    1.000000          1         1.000000             1
17   17     1.00000    1.000000          1         1.000000             1
18   18     1.00000    1.000000          1         1.000000             1
19   19     1.00000    0.990598          1         0.985013             1
20   20     1.00000    1.000000          1         1.000000             1
21   21     1.00000    1.000000          1         1.000000             1
22   22     1.00000    1.000000          1         1.000000             1
23   23     1.00000    1.000000          1         1.000000             1
24   24     1.00000    1.000000          1         0.878577             1
25   25     1.00000    1.000000          1         1.000000             1
26   26     1.00000    1.000000          1         1.000000             1
27   27     1.00000    1.000000          1         1.000000             1
28   28     1.00000    1.000000          1         1.000000             1
29   29     1.00000    1.000000          1         1.000000             1
30   30     1.00000    1.000000          1         1.000000             1
31   31     1.00000    1.000000          1         1.000000             1
32   32     1.00000    1.000000          1         1.000000             1
33   33     1.00000    1.000000          1         1.000000             1
34   34     1.00000    1.000000          1         1.000000             1
35   35     1.00000    1.000000          1         1.000000             1
36   36     1.00000    1.000000          1         1.000000             1
37   37     1.00000    1.000000          1         1.000000             1
38   38     1.00000    1.000000          1         1.000000             1
39   39     1.00000    1.000000          1         1.000000             1
40   40     1.00000    1.000000          1         1.000000             1
41   41     1.00000    1.000000          1         1.000000             1
42   42     1.00000    1.000000          1         1.000000             1
43   43     1.00000    1.000000          1         1.000000             1
44   44     1.00000    1.000000          1         1.000000             1
45   45     1.00000    1.000000          1         1.000000             1
46   46     1.00000    1.000000          1         1.000000             1
47   47     1.00000    1.000000          1         1.000000             1
48   48     1.00000    1.000000          1         1.000000             1
49   49     1.00000    1.000000          1         1.000000             1
> 
> save(muljj, file="jj.RData")
> proc.time()
[1] 1076.73   13.39 1129.70    0.00    0.00
> 
